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Abstract 

We simulate magnetic reconnection in the GEM problem using 
a two-fluid model with 10 moments for the electron fluid as well 
as the proton fluid. We show that use of 10 moments for electrons 
gives good qualitative agreement with the the electron pressure 
tensor components in published kinetic simulations. 

1 Overview 

This study is motivated by the following question: What is the sim- 
plest fluid model that can accurately replicate kinetic (PIC 
and Vlasov) simulations of fast magnetic reconnection in two- 
species collisionless plasma? Background for this question follows. 

A plasma is a gas of charged particles interacting with an electro- 
magnetic field. Simulations of plasma use a variety of models of plasma, 
which vary greatly in computational expense. For a given problem one 
seeks the computationally cheapest model that captures the phenom- 
ena of interest. We are specifically interested in the phenomenon of fast 
magnetic reconnection in "collisionless" (i.e., low-collision) plasma. 

We discuss the following sequence of plasma models for two-species 
(e.g. electron/proton) plasmas: kinetic: Vlasov/Boltzmann or PIC (particle- 
in-cell), ten- moment two- fluid plasma, five-moment two-fluid plasma, 
and MHD (magnetohydrodynamics, a one-fluid model). Each model in 
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this sequence can be regarded as a simplifying appoximation of its pre- 
decessor. The plasma community has used particle-in-cell (PIC) simula- 
tions as its standard first-principles plasma model. It is an approxima- 
tion of the Vlasov/Boltzmann model, which we take as the "truth". The 
Vlasov model is the collisionless version of the Boltzmann model. Colli- 
sionless versions of each model are hyperbolic and conserve entropy for 
smooth solutions, whereas collisional versions are diffusive and produce 
entropy. 

The Geospace Environmental Modeling (GEM) Magnetic Reconnec- 
tion Challenge is a benchmark magnetic reconnection problem. It was 
formulated to test the ability of plasma models to resolve fast magnetic 
reconnection in collisionless plasma. The initial GEM studies showed 
that PIC simulations exhibit fast reconnection, whereas reconnection is 
much slower in MHD (with simple resistivity) pQ. 

2 Plasma Models 

We implemented five-moment and ten-moment two-fluid models and 
studied their ability to match published results of PIC and Vlasov sim- 
ulations. Two- fluid models use gas-dynamics for each charged species. 
These fluids are coupled to Maxwell's equations by source terms. Hakim, 
Loverich, and Shumlak simulated the GEM problem with a five-moment 
two-fluid model [4] and Hakim simulated the GEM problem with a 
two- fluid model using 10 moments for ions and 5 moments for electrons 
[5], but we are unaware of any previous studies that simulate the GEM 
problem with a 10-moment electron fluid. 

Boltzmann/ Vlasov model. The Boltzmann equation asserts con- 
servation of particle number density / s (x, v,t) in phase space: 

ft/ s + V x -(v/ s ) + V,-(a/ s ) = C s ; 

here v = 7V « v is (proper) velocity, (where 7 = y/l + (v/c) 2 ~ 1 is the 
Lorentz factor), a = — (E + vxB) is the acceleration due to the electric 
field E and the magnetic field B, and C s is a collision operator which 
operates on the function (v,p) i->> / p (t,x, v), where p ranges over all 
species. The Vlasov equation (collisionless Boltzmann equation) asserts 
that C s = 0. The relations J = J2 S q s J v / s v and a = ^ s q s J v / s couple 
the Boltzmann equation to Maxwell's equations 

ftB = -VxE, V-B = 0, 

ftE = c 2 Vx5- J/e , V-E = a/e . 

Five-moment model. Generic physical equations for the gas-dynamic 
portion of the five-moment two-fluid model consist of conservation of 
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mass and balance of momentum and energy for each species: 
dtps + V-OsU s ) = 0, 

<9t(> s u s ) + V-(> s u s u s ) + Vp s = — p s (E + u s x B) + R s + V-a , 

d t £ s + V-(u s (£ s + p B )) + Vq s = V-(^-u) + J s E + u s R s + Q{ + Q*, 

where s is the species index (i for ions, i for electrons), p s is mass den- 
sity, u s is fluid velocity, £ s is gas-dynamic energy, q s = ±e is particle 
charge, m s is particle mass, and the species current is J s := (q s /m s )p s u s . 
For the pressure we assumed an isotropic monatomic gas: (3/2)p s = 
(£ s — p s ^g/2) . A linear isotropic entropy-respecting viscous stress clo- 
sure is a g = 2/i s (Sym(Vu s ) — V-u s I/3), where Sym denotes the sym- 
metric part of its argument tensor, I is the identity tensor, and p s is the 
shear viscosity. In these five- moment simulations, however, we neglect 
all collisional effects. So we neglect viscosity (a g = 0), heat flux (q s = 0), 
resistive drag force (R s = 0), resistive heating (Q( = 0), and interspecies 
thermal equilibration (Q* = 0). To couple these equations to Maxwell's 
equations we use the relations 

J = J^+J e , cr = (Ji + (j e , 

where a s = (q s /m s )p s is the charge density of each species. 

Ten- moment model. Generic physical equations for the gas-dynamic 
portion of the ten-moment two-fluid model consist of conservation of 
mass and balance of momentum and energy tensor for each species: 

dtps + V-(p s u s ) = 0, 

d t (p s u s ) + V- (p s u s u s + P 8 ) = — p s (E + u s x B) + R s , 
d t E s + 3V- Sym (u s E s ) - 2V- (p s u s u s u s ) + V-q s 

= ^-2 Sym (p s u s E + E s x B) + R s + 2 Sym (u s R s ) + Q* + Q*, 

where E s := J v / s vv is the energy tensor and P s := E s — p s u s u s is the 
pressure tensor. A linear isotropic entropy-respecting isotropization clo- 
sure is R s = (|(trP s )I — P s ) , where tr denotes the trace of its argu- 
ment tensor and for the isotropization period we used r s = tow det 5 Ps , 
which attempts to generalize the Braginskii closure; for the GEM prob- 
lem this means that Ti/r e ~ (mi/m e ) 5 ^ 4 . We neglect all other collisional 
terms: the heat flux tensors q s , the resistive drag forces R s , the frictional 
heating tensors Q(, and the temperature equilibration tensors Q*. 

For small viscosity (or fast isotropization) the viscosity is related to 
the isotropization period by p s ~ p s r s- We set To = 50. The ten- moment 
model offers the advantage of hyperbolic viscosity — that is, viscosity can 
be implemented without numerically expensive diffusive terms. 
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3 GEM Problem 



The GEM problem is posed on a rectangular domain with periodic 
boundary conditions in the horizontal direction and with conducting 
wall boundary conditions for the upper and lower boundaries. The ini- 
tial conditions are a Harris sheet equilibrium perturbed by "pinching" 
to form an X-point. 

Nondimensionalization. The GEM problem nondimensionalizes 
time by the ion gyrofrequency Qi = and the space scale by the ion 
inertial length ^, which is the distance traveled in an ion gyroperiod 
1/Vti at the ion Alfven speed va,i •= /iQ ^° no (where /io := (c 2 eo) _1 
is magnetic permeability). Under this nondimensionalization the model 
equations above remain unchanged with the exception that eo is replaced 
by e := 1/c 2 (where c is now the speed of light divided by the ion Alfven 
speed). 

Model Parameters. The GEM problem specifies that the ion/electron 
mass ratio is rrii/m e — 25 and the initial temperature ratio is Ti/T e = 
^Jrrii/m e = 5. The speed of light is not specified; we used the com- 
monly used value of 20, chosen to be sufficiently high to exceed other 
wave speeds but computationally feasible. 

Computational domain. The computational domain is the rect- 
angular domain [—L x /2, L x /2] x [—L y /2,L y /2], where L x = Stt and 
L y = 4-7T. The problem is symmetric under reflection across either the 
horizontal or vertical axis. 

Boundary conditions. The domain is periodic in the x-axis. The 
boundaries perpendicular to the 2/- axis are thermally insulating conduct- 
ing wall boundaries. A conducting wall boundary is a solid wall boundary 
(with slip boundary conditions in the case of ideal plasma) for the fluid 
variables, and the electric field at the boundary has no component par- 
allel to the boundary. We also assume that magnetic field runs parallel 
to and so does not penetrate the boundary (this follows from Ohm's law 
of ideal MHD, but we assume it holds generally). So at the conducting 
wall boundaries 




d y B x — By — d y B z — 0, 

E x = OyEy = E z = 0. 



Initial conditions. The initial conditions are a perturbed Harris 
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sheet equilibrium. The unperturbed equilibrium is given by 



B(y) = B tanh(y/X)e x , p(y) = —n(y), 



ni{y) = n e (y) = n (l/5 + sech 2 (j//A)), p e (y) = 



■p(y) 




■p(y)- 



On top of this the magnetic field is perturbed by 



SB = —e z x V(^), where 



ip(x, y) = ip cos(2ttx/L x ) cos(7ry/L y ). 



In the GEM problem the initial condition constants are 



A = 0.5, 



n = 1 



^0 = Bo/10. 



4 Method 



To simulate the ten-moment and five-moment systems we implemented 
a Runge-Kutta Discontinuous Galerkin solver with third-order accuracy 
in space and time on a Cartesian mesh. 

To suppress oscillations, after each time stage we limited the solu- 
tion in the characteristic variables of the cell average using a modifi- 
cation of Krivodonova's method (beginning with the coefficients of the 
highest-order Legendre basis polynomials and descending to lower order 
if limiting occurs) [8] . 

To clean the magnetic field we used a correction potential i/j in 
Maxwell's equations, as suggested in [7]: 

S £ B + VxE+xV^ = 0, ftE-c 2 VxB = -J/e , <9^+Xc 2 V-B = 0. 

These equations imply a wave equation that propagates the divergence 
constraint error V-B at the speed c%. We used x = 1.05. 

Since the GEM problem is symmetric, we imposed symmetry and 
solved the equations on the quarter domain [0,L x /2] x [0,L y /2]. 

5 Results 

We simulated the GEM problem with ten-moment and five-moment 
models and compared the results with the Vlasov simulations of [3] 
and the PIC simulations of [2 Q Their plots were made at a point in 

1 We have to negate some quantities because we call the vertical axis y and the 
out-of-plane axis z, opposite to the convention of [21 El - 
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Figure 5.1: Ramp-up of reconnected flux for 10-moment and 5-moment 
plasma. Fluxes are through the boundaries of the first quadrant of the 
domain. In each model the bottom three lines should coincide. Times at 
which the magnetic flux through the positive x-axis reached one nondi- 
mensionalized unit per model were 17.2 (our 10-moment) 13.0 (our 5- 
moment), 17.7 (Vlasov [3 ), 15.7 (PIC 0), 17.6 (10/5-moment [5 ), 15.6 
(5-moment [4]) and 15.3 (5-moment [6]). 



time when the flux through the positive x-axis approximately reaches 
1 nondimensionalized unitj^] This is shortly before the time when the 
reconnect ion rate peaks. 

As measured against kinetic simulations the ten-moment model re- 
connects at about the correct rate and the five-moment model reconnects 
a bit too quickly, perhaps because in the five-moment model compression 
in the outflow direction automatically causes increased pressure in the 
perpendicular directions, opening up the outflow region and artificially 
increasing the rate of reconnection. 

In contrast to the five-moment model the ten-moment model is ca- 
pable of representing an anisotropic pressure tensor. Our plots of the 
ten-moment electron pressure tensor components look like somewhat 
smudged versions of the corresponding plots for the Vlasov model. 
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Figure 5.2: Magnetic field when one unit of flux passes through the 
positive x-axis. These results compare well with Plate 1 of [2 and the 
magnetic field shown in FIG. 2 of [3]. 
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Figure 5.3: Diagonal components of the electron pressure tensor for 
10-moment simulation at flit = 18. These look like smudged versions of 
the plots in FIG. 5 of [3]. 
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Figure 5.4: Off-diagonal components of the electron pressure tensor for 
10-moment simulation at flit = 18. These look like smudged versions of 
the plots in FIG. 6 of [3]. 



less magnetic reconnection. Physics of Plasmas, 13, 092309 (2006). 

[4] A. Hakim, J. Loverich, and U. Shumlak. A high-resolution wave 
propagation scheme for ideal two-fluid plasma equations. J. Comp. 
Phys., 219:418-442 (2006). 

[5] A.H. Hakim. Extended MHD modelling with the ten-moment equa- 
tions. J. Fusion Energy, 27:36-43 (2008). 

[6] J. Loverich, A. Hakim, U. Shumlak. A Discontinuous 
Galerkin Method for Ideal Two- Fluid Plasma Equa- 
tions. Submitted to Physics of Plasmas, posted at 

\protect\vrule widthOpt\protect\href {http : //arxiv . org/abs/1003 . 4542}{arXiv : 1003 . 4542}vl 
(2010). 

[7] A. Dedner, F. Kemm, D. Kroner, C.-D. Munz, T. Schnitzer, and M. 
Wesenberg. Hyperbolic divergence cleaning for the MHD equations. 
J. Comp. Phys., 175:645-673 (2002). 

[8] L. Krivodonova. Limiters for high-order discontinuous Galerkin 
methods. J. Comp. Phys., 226:879-896 (2007). 



